Ginzburg-Landau Theory and Classical Critical Phenomena of Mott Transition 
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Theory of classical critical phenomena of Mott transition is developed for the dimensionality 
d < oo. Reconsidering a cluster dynamical mean-field theory (DMFT), Ginzburg-Landau free 
energy is derived in terms of hybridization function for a cluster-impurity model. Its expansion 
around a cluster DMFT solution reduces to a cj)'^ model. Inherent thermal Mott transition without 
spontaneous symmetry breaking is described by a scalar field reflecting the charge degrees of freedom. 
In the space of local Coulomb repulsion, chemical potential and temperature, a first-order transition 
surface terminates at a critical end curve. The criticality belongs to the Ising universality as a liquid- 
gas transition. Various quantities including double occupancy, electron filling and entropy show 
diverging responses at the criticality and discontinuities at the first-order transition. Particularly, 
electron effective mass shows a critical divergence in 2 < d < 4. Only at a certain curve on the 
surface, a filling-control transition and its related singularities disappear. We discuss detailed critical 
behaviors, effects of interplay with other critical fluctuations, and relevant experimental results. 

PACS numbers; 71.30.-l-h, 05.70.Fh, 05.70.Jk 71.10.Fd 



I. INTRODUCTION 

Spontaneous symmetry breaking occurs in a variety 
of phase transitions in correlated electron systems. In 
particular, at partial but commensurate electron fillings 
n, some of them appear as a metal-insulator transition 
(MIT). For an integer n, the effect is most pronounced, 
because the MIT in this case is triggered by severe com- 
petition between kinetic energy and local Coulomb re- 
pulsion U for localized orbitals. This MIT is classified 
into a prototype of Mott transition^. Here, translational 
invariance is partly broken and a unit cell is enlarged. 
This results in a reconstruction of the bands which are 
either fully occupied or empty in the reduced first Bril- 
louin zone, and hence an insulator explained by Slater^. 

Actually, Mott transition takes place with or without 
symmetry breaking. It was postulated as a phase tran- 
sition that occurs at partial but integer n from metal to 
(magnetically) long-range ordered or disordered insulator 
with local spin momentsi. A central issue in this paper is 
the latter prototype of Mott transition occurring without 
any spontaneous symmetry breaking, which is henceforth 
referred to as inherent Mott transition; a first-order Mott 
transition curve separating the metallic and the insulat- 
ing states terminates at a critical end point, leaving a 
crossover at higher temperatures, as shown in Fig. ^ (a-) • 
This insulating state is realized by a Hubbard-band split- 
ting due to a large value of U compared with a bandwidth 
W^. The first-order nature may be assigned to elec- 
tron correlations even without lattice degrees of freedom. 
An analogy of the inherent Mott transition to a liquid- 
gas phase transition^ was first discussed by Cyrot^ and 
Castellani el al& and demonstrated by a number of exper- 
iments on d-electron systems such as V2O3 and NiS2iii 
as well as molecular materials /t-(BEDT-TTF)2)<S. 

In the proximity of Mott transition, there appear im- 
portant phenomena like high-Tc cuprate^i or organic su- 



perconductivity (SC) and possibly spin liquid^" in the 
resonating valence bond state-'^^, in addition to antifer- 
romagnetism (AF) due to the superexchange mechanism 
within a Mott insulating state. Figures ^ (b-d) show 
schematic phase diagrams for interplay of the inherent 
Mott transition with other phase transitions, which will 
also be analyzed in this paper. In (b), the inherent 
Mott transition is replaced by the Slater's mechanism 
with a symmetry-breaking Mott transition which has 
a larger scale of temperature T than the Mott criti- 
cal temperature T^. In (c) and (d), below Td-, there 
occurs interplay of the inherent Mott transition with 
other symmetry-breaking transitions that do not appear 
as an MIT. These phase diagrams are relevant to (c) 
NiS2-a;Se:c and i?Ni03 (R: La to Lu)i*i for interplay with 
AF and «;-(BEDT-TTF)2X for interplay with SC and 
AFi^, (dl) Va-xCr^O^ii and k-(BEDT-TTF)2X (under 
a magnetic field)i^ for interplay with AF, and (d2) k- 
(BEDT-TTF)2Cu2(CN)3ia or i?2Mo207 {R: Nd to Y)il 
for interplay with SC or ferromagnetism. 

On the paramagnetic metallic side of the Mott tran- 
sition, quasiparticle mass and hence a specific heat co- 
efficient C/T are enhanced and even diverge at quan- 
tum Mott criticality^, as obtained in infinite dimensions 



d 



The quantum Mott transition for the single- 



band Hubbard model was described as a holon-doublon 
binding- unbinding transition"'^^. Holons and doublons are 
bound by a large value of U /W at half filling, and then 
no longer propagate freely, yielding a Mott insulating 
state. The binding energy gives a Mott gap and thus 
the inverse of the localization length ^loc. This ^loc cor- 
responds to the holon-holon (or doublon-doublon) and 
the holon-doublon coherence lengths. However, under- 
standing critical properties of Mott transitions remains 
open in realistic spatial dimensions, whether it is quan- 
tum or classical. Now, it is important to clarify their 
universal properties. 
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FIG. 1: Schematic phase diagrams for (a) inherent (liquid- 
gas) thermal Mott transition and a related crossover, and (b- 
d) the interplay with other symmetry-breaking phase tran- 
sitions, (b) The inherent Mott transition is masked by the 
Slater-type one. (cl/c2) There occur two additional phases 
separated by the first-order transition from each other, which 
extends over as an inherent Mott transition at high temper- 
atures. (dl/d2) In the metallic/insulating side, there occurs 
a symmetry-breaking phase transition. Fluctuation-induced 
first-order transitions possibly appear in the vicinity of meet- 
ing points of these curves-. 



Second-order and some first-order classical phase tran- 
sitions can be understood by analyzing Ginzburg-Landau 
(GL) free energy. We can even describe critical phenom- 
ena of conventional second-order ones in terms of univer- 
sality class uniquely specified by the number of compo- 
nents of relevant order parameters in each system dimen- 
sion, as far as the interaction is short-rangei-:. Then, it 
is fundamental to find a minimal complete set of relevant 
order parameters. For symmetry-breaking phase transi- 
tions, one may be ready to select the order parameters 
conjugate to the symmetry-breaking fields and hence to 
construct the GL free energy. In contrast, for liquid- 
gas phase transitions occurring without any symmetry 
breaking, the two states on the both sides of the first- 
order phase transition are adiabatically connected. The 
relevant order parameter may be hidden, and hence it is 
nontrivial to derive the GL theory. 

Only in c? = oo, GL theory of Mott transition has been 
developed by Kotliar el ali^ along the formalism of the 
dynamical mean-field theory (DMFT)iSiS£. "Order pa- 
rameter" is derived from hybridization function in the 



single-impurity Anderson model as a bosonic scalar field. 
Its finite-temperature criticality belongs to the mean- 
field (MP) universality. This GL theory together with 
DMFT solutions of the Hubbard model implies that the 
doublon and the charge susceptibilities 

1 d^F 

- 

simultaneously diverge at Mott critical end curve in the 
space of {U, fi, T)^ where F, fjL and N are a free energy for 
the grand canonical ensemble, a chemical potential and 
the number of atomic sites, respectively. However, de- 
tailed analyses in the whole {U, /i, T) space remain open. 
Besides, in finite dimensions d < oo, even derivation of 
GL free energy is missing. 

In this paper, we derive GL free energy of Mott tran- 
sition for d < oo. GL free energy describing the inherent 
thermal Mott transition turns out to be similar to that 
of the Ising model under a magnetic field. Then, we con- 
firm that the analogy to the liquid-gas phase transition^ 
holds, but with crucial modifications of the DMFT re- 
sults by spatial correlations^ in 2 < d < where d,, = 4 
is the critical dimension^^. This is consistent with the 
absence of thermal Mott transition for d = 1 , a phase di- 
agram proposed for the Hubbard model in d = 324, and 
a recent experimental scaling analysis of the resistivity 
in Vs-.Cr.OgS^. 

We propose a {U, fi, T) phase diagram for the inherent 
Mott transition (Fig. EJ. A first-order Mott transition 
surface (shaded area) terminates at a critical end curve 
(thick solid curve). Both bandwidth-control and filling- 
control Mott transitions occur there, except at a special 
curve (thick dotted curve) and its end point ([/„•, Acn ^cr) 
(open circle) where a bandwidth-control transition takes 
place but not a filling-control with U being fixed. The 
double occupancy D is discontinuous at the first-order 
transition surface and xdd diverges at the critical end 
curve. The electron filling n also shows a discontinuity 
at the surface but at the special bandwidth-control curve 
where its singular component vanishes. The same critical 
divergence as xdd emerges for x„„ and 

1 d^F 

at the critical end curve but at (C/cr, Acr, Tcr) where di- 
verging components of Xnn and Xon vanish. Xon oo 
at criticality is a quite new feature. These refiect that 
each of D and n contains a linear term in the relevant 
scalar field. Variations of C/, /i or T merely change the 
expansion coefficients of the GL free energy. Taking into 
account this aspect, we have succeeded in obtaining a 
unified picture of the Mott transitions with respect to 
different control parameters. 

In Sec.m to establish a basis of GL arguments for Mott 
transition, we reconsider a general theoretical framework 
of a cluster extension of the DMFT. 
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FIG. 2: Schematic phase diagram of the inherent Mott tran- 
sition in the {U,fi,T) space. A first-order transition surface 
(shaded surface) terminates at a critical end curve (thick solid 
curve). Within the surface, there exists a curve (open dotted 
curve) where a first-order bandwidth-control transition with 
a discontinuity of D occurs but not a filling-control. This 
special curve terminates at a critical end point (C/cr, Acr, Tcr) 
(open circle) on the critical end curve. At this point, a second- 
order bandwidth-control transition with Xdd — > oo occurs. 
Except at the special curve and point, both D and n show 
discontinuities at the first-order transition surface and their 
U and ^ derivatives diverge_ at the critical end curve. The 
dashed line represents U — Uci and T — Tcr, while the dash- 
dotted line /i = /icr and T — Tcr- 



In Sec. mil GL free energy functional for this cluster 
DMFT is constructed in terms of hybridization function 
of a cluster-impurity model. In the limit of an infinite 
cluster size, expanding the free energy around a critical 
point and transforming fermionic field to bosonic, we ob- 
tain a (f/^ model for the GL free energy. 

In Sec. IIVI we consider the inherent Mott transition 
with finite-temperature criticality. We reveal universal 
properties of various physical quantities, the connection 
among Mott transitions driven by different control pa- 
rameters, and the phase diagram shown in Fig. |21 

Interplay of the inherent Mott transition with other 
symmetry-breaking transitions is a topic of Sec.0 From 
GL arguments, as far as critical regions of the other tran- 
sitions do not overlap, the Mott criticality appears at an 
end point of the liquid-gas phase transition or is masked 
by the other phase transition, as shown in Fig. ^ On 
the other hand, one-loop renormalization group analyses 
show that in the presence of the appreciable overlap, it 
is replaced with a fluctuation-induced first-order transi- 
tion, as in a previous discussion on (c) for interplay with 
SO and AFi^. 

Lastly, we summarize the present results and comment 
on problems left for future studies in Sec. IVII 



II. CLUSTER DYNAMICAL MEAN-FIELD 
THEORY 

There have been attempts to restore spatial correla- 
tions ignored in the OMFTiSlSLSi. On one hand, cor- 
relator projection methodSi combines the DMFT with 
projection methods^^'^^-^^ of electronic operators leading 
to continued-fraction expansion of Green's functions in 
energy uj^. On the other hand, cluster approaches^&Si 
solve a cluster problem containing multiple impurities in- 
stead of a single impurity in the DMFT. They should be- 
come exact if the cluster size is taken infinite. However, 
depending on the formalisms. Green's functions and self- 
energy parts have different properties. In this section, 
we reexamine a general framework of the cluster DMFT, 
which gives a basis for constructing the GL free energy. 



A. General requirements 

We impose the following three requirements, whichever 
the cluster is introduced in a real or a momentum space. 

(i) Dynamical quantities such as Green's function and 
self-energy parts are manifestly causal. Difficulty in 
satisfying this has already been overcome in the liter- 
atureSS^. 

(ii) Each symmetry of the Hamiltonian is sponta- 
neously broken in a solution if and only if order parameter 
conjugate to the symmetry breaking field remains finite, 
except the translational invariance which is not necessar- 
ily guaranteed because each cluster is treated separately 
from the others. This requirement is crucial for classify- 
ing phase transitions. 

(iii) Self-energy parts including the secular MF 
(Hartree-Fock and/or BCS) terms are only taken into 
account within the cluster. This facilitates arguments 
to derive a cluster DMFT equation. Although this re- 
quirement introduces an unphysical barrier between the 
cluster and the exterior medium, this boundary effect is 
0(Lf~^) for the cluster size Nc = and hence negligible 
compared to as Lc ^ oo at least for d> 2. 



B. Lattice model 

To be explicit, we consider an Af -component electron 
system described by the following Hamiltonian; 

= £cri,o-2(k)c^^^Ck,cr2 

k,(Ti ,(T2 

+ Y Y ^'^i:'^2;'T3,<T4(ki,k2;k3,k4) 

ki ,■ ■ ■ ,k4 (71 ,-■ ■ ,(74 

^''ki,(7i'^k2,(72''''3>'^3'"l<^<l''^<l' 

Here, Ck,<7 and cj^ ^ are annihilation and creation oper- 
ator of electron with the momentum k and the internal 
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component degrees of freedom a. For instance, a rep- 
resents spins, orbitals and sublattices. V represents a 
short-range Coulomb interaction for the electrons and is 
finite only when ki + k2 = + k4. E^^^a^ (k) is given by 



,(k) = -t^i^<^2(k) - ^(5^1, 



(5) 



tcri,(T2(k) is an electron transfer with the momentum k 
from CT2 to (Ti and satisfies X)k ^o'i,o-2 (k) = 0. denotes 
a chemical potential, and /icti,<t2 represents a traceless 
uniform field which is locally coupled to electrons and 
breaks a symmetry in the component degrees of freedom, 
e.g., magnetic field, energy- level splitting among different 
orbitals, lattice distortion in the Jahn- Teller term, etc. 

We define electron thermal Green's function at inverse 
temperature /? = 1/T as 



/•/3 

a,,.(ic^„;k,k') = - dT ([ck,.(r),4, J+)e'^-^ 
Jo 



(6) 



with [A,B]+ = AB + BA and Ck,a(r) = e'^^Ck.ae"^^. 
Henceforth, uj„ ~ {2n + 1)'kT and flm = 2m7rT de- 
note fermionic and bosonic Matsubara frequencies, re- 
spectively. The Dyson equation takes a matrix form; 

G'~"^(ia;„; k, k') = iw„/ - e(k) (5k,k' - S(iw„; k, k'). 

We have introduced notations G, E, / and e for M x 
M matrices Ga-,cr', ^a-.a-', 5a.<j' and e^.o-", respectively, 
where Ecr.cr' k, k') is the self-energy part. 



C. Cluster model 

To calculate S(«w„; k, k'), we introduce as a represen- 
tative of the original problem a cluster model consisting 
of Nc sites in the real space as in the cellular DMFT^^ 
or Nc momentum cells as in the dynamical cluster ap- 
proximation (DCA)^. We define the Green's function 
Gc{ii^n] i, i') and the self-energy part Sc(iti;„; «, i') for the 
cluster model, on which we impose self-consistent condi- 
tions between cluster and lattice variables; 

S(zc^„;k,k') = -^^i?(k,^)^],(^w„;^,^')i^^(^',k') 



(8a) 



GS^^-i,i') = l^i?t(z,k)G(za;„;k,k')-R(k',z'). 

k,k' 

(8b) 

Here, i and i' represent the Nc degrees of freedom of 
the cluster, i? is a transformation matrix from cluster to 
lattice variables characterizing the way of introducing a 
cluster. It obeys the orthonormal condition. 



l^i^t(^,k)i^(k,^') 



(9) 



If we take i?(k, i) = for k belonging to zth mo- 

mentum hypercube centered at k(i) and otherwize, 
it gives a modification of the DC A^^i?^ . The choice of 
i?(k, i) = e*'"'(') with r(i) being the position of ith site 
within the cluster is taken in the cellular DMFT^^. 

The cluster model can be derived by integrating out 
the remaining degrees of freedom other than the cluster 
along the concept of the cavity construction^*^. However, 
in finite dimensions d < oo, the exterior degrees of free- 
dom mediate additional one-body and many-body inter- 
actions among the cluster degrees of freedom^SiSl. Ne- 
glecting the exterior- mediated many-body interactions, 
we only include the one-body term arising from the hy- 
bridization between the cluster and the exterior as in the 
literature^SiSl. 

Thus, we obtain the following path-integral expressions 
of the partition function and the effective action for the 
cluster in terms of Grassman field / and / representing 
the electronic degrees of freedom within the cluster. 



Zc = JvfVfexp[~Sc 



(10a) 



Sc ^ -f3 ^ /j^,<^,(itj„)©^/_^^(itJ„,ii,i2)/j2,CT2(«^«) 



n,(Tl ,(72 ,«1 ,22 

dr Va^ ,a2;a3,a4{il, 12; is, U) 







^ fil ,cn ir)! 12, <y2 (^)/i3,<''3 ij)fii,ai (j) 



(10b) 



with 



,(^2 ;C3 ,'74 

= ^ ^o-i,<J2;(j3,<T4(ki,k2;k3,k4) 

Xi?t(ki, ii)i?t(k2, I2)i?(k3, h)R{ki, ii). (11) 

© plays a role of the Weiss field acting on the cluster in 
the original lattice model, and is now nonlocal within the 
cluster in contrast to the conventional DMFT where it is 
completely local. © is obtained from 

©"^(icj„; i, i') = iujn5i,i'i - £o(«, «') ~ A(icj„; i, z'), (12) 

where 

eo(«,0 = ^5]i?t(z,k)£(k)i?(k,z'). (13) 

k 

is the intra-cluster part of £(k), and 

A(itj„;i,i') = -^^£|(i,k)Ge(«t^n;k, k')ei(k',i') 

k,k' 

(14) 

is the exterior-mediated dynamical one-body interaction, 
i.e., the hybridization function in the context of a multi- 
impurity extension of the Anderson impurity model. 
Here, we have introduced the Green's function for the 
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original lattice model in the absence of the cluster 
G'c(«w„;k, k') = G'(icj„;k, k') 

-772 E E E^(*^";i^'''i)^(''i'*)^^(*''^2) 

ki ,k2 k'j jkj 

xG'-i(ic^„; k2, k^)ii'(k'2, k'i)G(iw„; k'^, k'), 

(15) 

and hybridizations between the cluster and the exterior 
elii,^) = i?t(i,k)£(k)-^eo(i,i')^^(*',k),(16a) 

ei{k,i) = e(k)i?(k,i) - ^i?(k,i')eo(i', «)• (16b) 

The second terms in Eqs. (|16a|l and Ijl6b|l remove the 
intra-cluster contributions to the hybridizations to avoid 
a double counting of the cluster degrees of freedom. 

The Green's function Gc and the self-energy part Ec 
for the cluster are obtained by solving the cluster model 
Eq. 110|l and are related through the Dyson equation 

G^^(it^„;i,z') = (S"^(«a;„;i,i') - Sc(icj„; i, i')- (17) 

Finally, Eq. Ij8b|l is reproduced by substituting Eq. I|14j) 
together with Eqs. ifT^ and (|Sa|l into Eq. ifT^ . and then 
the result into Eq. (|17|l . Equation together with 
Eqs. (jHl and ^ defines a set of cluster DMFT self- 
consistent equations. 



III. GINZBURG-LANDAU FREE ENERGY 
FUNCTIONAL 

In this section, we construct a GL free energy func- 
tional in terms of the hybridization function so that for 
symmetry-unbroken states on the Bethe lattice in d = cx3, 
it reproduces the Kotliar's result In d = 00, the hy- 
bridization function was regarded as a metallic order pa- 
rameter. For d < 00, it is nonlocal within the cluster. 



Therefore, the cluster-average low-energy behavior of the 
trace of the hybridization function plays a role of the 
metallic order parameter. Then, we define a bosonic field 
corresponding to the "order parameter" and describing 
the soft mode at criticality, and finally obtain a 0* model, 
which facilitates analyses of the GL free energy as will be 
done in Sees. IIVI and Ivl 



A. Construction of free energy functional in terms 
of hybridization function 

As is clear from Eqs. (|1U|) and (|12|l . the effective action 
Sc and the partition function Zc are functionals of A in 
addition to a set of model parameters denoted by s = 
{U, t, t', fi,h,T, - ■ ■ }. Then, the cluster Green's function 
Gc and self-energy part Ec sue determined by A from 

- ~ 1 5F - 

Gc[A,s\{iuj„-i,i') = --f [A,s](^c^„;^,^') (18) 
T SA 



F,[A,s\ = -TlogZ,[A,s], 
and from Eqs. H12|l and (|17|l . namely, 



(19) 



T.c[A,s\{iuJn]i,i') = iUnSij'I - eo{i,i') - A{iuj]i,i') 
-G-l[A,s1(^t^„;^,^'), (20) 

respectively. Therefore, the lattice Green's function 
Eq. (UI) with the self-energy part Eq. l(5al) can be regarded 
functional of A as well; 



G-MA,s-](iw„;k,k') = iujj-e{k) 



5k,k' 



J2 ^(k, »)Sc[A, s\iiu;„; i, i')R^{i', k'). (21) 



Finally, via a Legendre transformation from Gc to A, we 
can construct GL free energy i^cL as a functional of A; 



5Fgl[A,.§1 = (5Fc[A,s1-- Tri?t(,',k)G[A,s1(io;„;k,k')i?(k',z)5A(ic^„;i,z'), (22) 

riMM.' 

SO that the stationary condition directly gives the cluster DMFT equation given by Eq. ((SbJ; 

f [A, s](zc.„; z, i') ^ [A, s] = Gc[A, s\{ilo^; i, " ^ E ^^(*' ^)^[A, s\{iu;^- k, k')i?(k', i') ^ 0. (23) 

k,k' 

I 

Especially for the Bethe lattice in d = 00, Eq. 114|) i = Ymi^^^t/ ^/d and the nearest-neighbor transfer 
is simplified as A{iuin) = PGc{iuJn) with Nc — 1, and hence Eas. I|22|l and H23I) reproduce the Kotliar's re- 
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Fgl[A,s] =F,[A,sl-T^A(ic^„)V2P. (24) 



It is remarkable that A plays a role of an auxiliary field 
in the present theory. In conventional functional-integral 
techniques like the Hubbard-Stratonovich (HS) transfor- 
mation in the Hubbard model, auxiliary fields can be 
introduced for describing the magnetic order parameter 
and long-range ordering. There are three important dif- 
ferences between these two representations. 

(i) The saddle-point solution leads to a cluster DMFT 
in the present case, while it only gives the Slater's MF 
theory in the conventional functional-integral method. 

(ii) In the present theory, A includes the charge degrees 
of freedom and their instability in addition to the mag- 
netic ones, while random phase approximations (RPA) 
or (/)4 theory for itinerant magnetism based on a simple 
HS transformation^Si^SiSi can not properly describe the 
charge degrees of freedom when their fluctuations develop 
near the Mott criticality. 

(iii) In the present theory, the auxiliary field is 
fermionic and contains momentum-resolved fermionic 
spectral properties which gives crucial contributions to 
the momentum dependence of the self-energy correction 
in low dimensions, while in the above theories for itin- 
erant magnetism, one only deals with an bosonic field 
which is locally coupled to fermions. 



B. Expansion of Ginzburg-Landau free energy 
functional around a critical point 

In order to obtain GL expansion of the free energy, we 
assume that in the original model, a criticality appears at 
a set of model parameters s = Sa (temperature Tcr, lo- 
cal Coulomb interaction Ucr, chemical potential /icr, etc). 
We take the limit of an infinite cluster size A^c ^ c» so 
that Nc/N becomes a finite value less than unity. We fur- 
ther assume a periodic boundary condition. Then, the 
hybridization function A(iw„,z,z') can be expressed in 
the energy-momentum space (iw„,k) as A(iti;„,k), since 
due to the translational symmetry, the momentum con- 
servation law as well as the energy conservation holds. 
(When we consider a translational-symmetry-breaking 
transition with a unit-cell enlarging, we need to take this 
enlarged unit cell at the initial stage of Eqs. Q and l^.) 
Then, expanding the free energy _Fgl[A,s| around the 
cluster DMFT solution Acr at the critical point Scr in 
terms of 5 A = A — Acr, we obtain 

Fo.=Fi°l+Fi^l+Fi'l+Fil+F^t + --- (25) 



with a regular part Fq-^ = FGL[Acr, s\ and 

4l = T^Tr(f(i)[A„,s1(fc)(5A(fc)) (26a) 

k 

^GL = 7 E r(^l|...i;.JA„,.in'5AK,-(26b) 

Ki.-.Ki 3 = 1 

Here kj and Kj are abbreviations of (iwn^. ,kj) and 
(aj ,(Jj)(^ {ioJuj , ) J respectively, and we have introduced 

rS,,...,..[A,.l.^(n^)rglA„^(27, 

for £ > 2. In principle, this enables us to obtain the 
GL expansion parameters for explicit models, although, 
practically, it may not be easy since it requires extrapo- 
lations of cluster Green's function, self-energy part and 
vertex parts into the limit Nc ^ oo with Nc/N < 1 being 
fixed. 

C. ^^-model representation 

We label the real eigenvalues of the hermitian 
matrix r^^^^^^,|^^^^, [A„., ^(fcijfca) ee r^j^^JAcr,s] as 

{7^^-'(p)[s]} with A — 0, • • • , — 1 and an abbrevia- 
tion p ~ (ir2„i,p), and the corresponding eigenbases as 
{'i'a^a'\x[s\ik\p)}. Namely, 

E rif,.;|..,.^[Acr,«1(fcl|fc2)*.„.i|A[.1(fc2b) 

fT2 ,(72 

^^,,y^\^mk^\pW^\s\ip). (28) 

We further introduce 0a (p) as the component of 
(5Ao-^cr'(fc) projected onto 'i'cr,a'\x[s\{k\p); 

SA,^,,{k) = ^^'.,.'|A(fcb)0A(p). (29) 

A,p 

Owing to the energy-momentum conservation, 
^o-,(T'|A[s](fc|p) contains an energy- momentum 
k — p — {iujn — irimjk — p). Namely, the fermionic 
function A^^o-' is expressed as a convolution of the 
fermionic function ^'^^o-'ia and the bosonic function 

Substituting Eq. ((23) into Eqs. (j^^ and (|26b|) . we 
obtain 

4l = r7i'VA(0) (30a) 

Fol = ^J2^^x\p)\MPr (30b) 
p 

^GL = f E tSps.Ps'^-Pi'^^^M (30c) 

^GL = ^ E ^plp2,P3,Pi'i^Pl'l'P2<t>P3<l>P4 (30d) 

Pu--,Pi 
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where Pi denotes an abbreviation for a set of Xi and pi 
and we have also introduced 



K 



(7,cr' I A 



mm, 



(31) 



TpI,. 



.Pi 



A. Ising universality class 

It is reasonable to assume that there occurs an inher- 
ent thermal Mott transition for c? > 2 that does not break 
any symmetry spontaneously, as in liquid-gas phase tran- 
sitions. Our classical GL free energy reads 



GL 



(32) T 



with £ ~ 3 and 4. The energy-momentum conserva- 
tion guarantees that 7'^^ has only a static uniform com- 
ponent and that ... vanishes unless X]j=iPj — 



(33) 

h, r and u represent the first-, the second- and the fourth- 
order coupling constants, respectively, which are func- 
tions of s. We have omitted the zeroth-order term F, 



(0) 
GL' 



I]j=i(if^mj,Pj) = for £ 

„(3) .^A .,(4) 



3 and 4. We note which is regular and does not alter singular behavior ob- 



that P2 P3 and 7pj^ p^ p^ are symmetric under any 
permutation between Pj and P^/ , and that 7^^^,' = 
W|m"|m"(^'i'P2,P3) = holds in the presence of 0{m) 
symmetry in with m' = mg, • • • , TOq + m — I and 
< m" < AP - 1. Eqs. (EOJl together with Eqs. lO and 
are reduced to a 0"* model. 
Apparently, phase volumes of the components (c, c') 
and the energy-momentum (iti;„,k) for (5A(ia;n,k) coin- 
cide with those of A and (il7m,p) for 0A(«f^m,p), re- 
spectively, because we have only performed a unitary 
transformation Ea. (|29|l . Among these degrees of free- 
dom, there exists a unique field (/<o(*^miP) whose static 
uniform part need not spontaneously break a symmetry 
of the Hamiltonian and can be directly related to the to- 
tal electron filling and the total multiple occupancy. In 
particular, if the Mott transition occurs without any sym- 
metry breaking in the single-band model, then the com- 
ponent degrees of freedom (cr, a') and A disappear in all of 
Green's functions, hybridization functions and the field 
(f)\{p)- Therefore, a field intrinsic to Mott transition has 
only one component, which is referred to as (/^(p) = 4>q{p)- 
This result is physically natural since this scalar field rep- 
resents nothing but the charge degrees of freedom. Other 
field components describe relative changes for the spin, 
orbital, intra-sublattice, and/or pairing degrees of free- 
dom. 



IV. CLASSICAL CRITICAL PHENOMENA OF 
INHERENT MOTT TRANSITION 

In this section, we consider the inherent thermal Mott 
transition in the absence of interplay with other degrees 
of freedom. The GL free energy in terms of the scalar 
field (f{p) yields the Ising universality for 2 < d < dc 
due to critical order parameter fluctuations and the MF 
universality for dc < d a.s in d = oo^^ where dc = 4 is 
the critical dimension. Here, we discuss various universal 
properties characteristic of the Mott criticality 2 < d < 
oo. We clarify a {U, fj,, T) phase diagram shown in Fig. 12 
and the connection between the bandwidth-control and 
the filling-control Mott transitions. Experimental tests 
of observable quantities are also proposed. 



("3") 

tained below. The third-order term has been re- 



moved by redefining (p{p) as 



7o AO (Pl ' P2 , P3)7o?0,0,0 (Pl , P2 , P3 , P) ■ 



„(4) 



Since the dynamics does not alter the universality of ther- 
mal transitions and the nonlocality of the quartic cou- 
pling u is irrelevant, they have been neglected in Eq. I|33|l . 

In this GL free energy, a criticality takes place at 
(/i, r) = (0, rcr(it)) where rcr(u) is a counter-term can- 
celing out the rcnormalization of u to r. Then, the criti- 
cality belongs to the Ising universalitjs^; 



S.\r,h=0 « |r - rcr(u)|"'', 

Ch\r,h=Q oc |r-rcr(u)r" 

(V')|rJi=o « {rcr~rf for r < rcr, 

x\r,h=o oc |r - rcr(M)r'^, 

cx sgn{h)\h\'l\ 



(34a) 
(34b) 
(34c) 
(34d) 
(34e) 
(34f) 



where ^ is the correlation length characteristic of the 
fluctuations of iy9, Cu is a speciflc heat at a constant 
h on the metal-insulator boundary = 0, and x = 
— j^d'^F/dh? is the susceptibility of ip. Scaling hypothe- 
sis gives relations among the exponents; a = 2 — dv, (3 = 
iy{d-2+r])/2, 7 = i^{2~T]), and 6 = {d+2-r]) / {d-2+r])^. 
In the case oi d = 3, ^ 0.64, a - 0.09, (3 - 0.33, 
7 ^ 1.25, and S ^ 4.8. In the case oi d — 2, ly = 1, a = 0, 
(3 — 1/8, ^ — 7/4, and S ^ 15. Furthermore, varying 
h across the Mott metal-insulator boundary h = gives 
a flrst-order transition for r < rcr{u) and a crossover 
for r > rcr(u). However, when we consider the scaling 
properties in the s space, the absence of the inversion 
symmetry if —ip bears complicated properties as in 
liquid-gas phase transitions^*^, which will be discussed 
later. 



B. Phase diagram 

Next, we discuss universal behaviors in terms of the 
original model parameters s. Mott metal-insulator 
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boundary can be defined as the trajectory of h[s\ = 0. 
The distance from the criticahty is determined by {h, 6r) 
with 5r[s\ = r[s\ ~ rcr{u[scr*[s\]), where Scr*[s| is such Sa 
point as gives a minimum value of r[s| — r(u[s*cr]) if Scr is 
not restricted to one point in the s space. This imphes 
that for instance, when s has three components, there 
appears a first-order transition surface terminating at a 
critical end curve. This is consistent with the recently 
proposed {U,t',T) phase diagram of the two-dimensional 
square-lattice Hubbard model24 where t' is the second- 
neighbor transfer. Therefore, we can fix the transfer pa- 
rameters and consider changes in U and T. We can also 
fix other short-range Coulomb interaction parameters. 

In addition to U and T, the chemical potential fj, plays 
a crucial role in Mott transition, since it controls the 
electron filling n. Quantum Mott transition should take 
place at a partial but integer filling, and hence there are 
two MIT boundary values of /i which correspond to the 
edges of the Mott gap. This boundary gives a curve in 
the (/i, U) space. At finite temperatures, a choice of n 
uniquely determines fi. Then, the Mott transition ex- 
tends to small electron and hole doping regions at T > 0, 
when it survives as a thermal phase transition. 

When the quantum Mott transition in the ([/, /i) space 
has a second-order character, then critical fiuctuations 
may immediately smear out the thermal transition in 
the space of {U,fi,T). On the contrary, when the quan- 
tum Mott transition has a first-order character, then 
it is robust against the thermal fluctuations and pro- 
duces a first-order Mott transition surface terminating 
at a finite-temperature critical end curve. This is be- 
cause in the ground states on both sides of the first- 
order transition, there exist neither critical fluctuations 
nor strong classical fluctuations characterized by power- 
law or exponential temperature dependence of the corre- 
lation length. Although such strong thermal fluctuations 
yield an abrupt collapse of long-range order containing a 
finite order parameter amplitude, these fluctuations ap- 
pear above the quantum second-order transition in the 
one-dimensional Ising model and in the two-dimensional 
Heisenberg mode]~, but not at the first-order transition 
in the present case. 

Concludingly, the finite-temperature inherent Mott 
transition is schematically described by a phase diagram 
shown in Fig. 12 Here, the first-order transition sur- 
face corresponds to a trajectory of h[U,fj,,T] = 0, while 
the critical end curve is further specified by r[U,fj,,T] = 
rcr{u[U, fJ.,T]). Similar first-order Mott transition surface 
terminating at a finite-temperature critical end curve has 
also been discussed in infinite dimensions'^, although we 
will find novel crucial aspects of critical thermodynamic 
properties in the rest of this section. 

To recognize the shape of the Mott transition boundary 
and to discuss critical properties in terms of s, it is impor- 
tant to note that h[s\, r[s\ and u[s\ are all regular func- 
tions, which is guaranteed in the GL arguments. There- 
fore, we can safely expand {h,r) around (0,rcr) in terms 
of s — s^j. For instance, dh/dfi is well defined around 



the criticality. Then, as depicted in Fig. |51 there exists a 
curve (open dotted curve) determined by h = dh/d^ = 
with r < Tc, where an infinitesimal increase of U but not 
any change of /i drives the system from metal to Mott in- 
sulator. Namely, there occurs a bandwidth-control Mott 
transition, but not a filling-control with U being fixed. 
The shape of the Mott transition surface at finite temper- 
atures is characterized by dU/dfj. = —{dh/d^)/{dh/dU), 
which vanishes toward the curve. This means that at 
least at T > 0, the surface is free from singularity and 
particularly rounded around the curved. Only at T = 0, 
there remains a possibility of a nonsingular behavior 
in the shape of the Mott phase boundary surface with 
respect to /i, reflecting the Fermi distribution function 
which appears in expressions for h, r and u. 

In most of the cases including d — oc^^ and d = 324, 
the first-order Mott transition surface is not parallel to 
the T axis in the {U, ^, T) space and a metallic state ap- 
pears at the lower temperatures. This is because stronger 
low-energy quasiparticle excitations in the metallic side 
contribute to a smaller entropy than in the insulating 
side. As fluctuations of other degrees of freedom like 
spins than charges develop, the entropy difference be- 
tween the metal and the insulator decreases and hence 
the surface becomes more vertical. 

More generally, locations of /i = and r = rcr(u) in 
the s space are nontrivial. For instance, they are not 
necessarily perpendicular to each other even around the 
criticality, in contrast to the usual Ising model under the 
magnetic field. Namely, it may be difficult to find the 
two routes of /i = and r = rcr(u). This may complicate 
detailed experimental scaling analyses, as in liquid-gas 
phase transition. 

C. First derivatives of free energy 

Here, we reveal singular behaviors of the double occu- 
pancy D, the electron filling n and the entropy per site 
S. Detailed expressions for general first derivatives of 
the free energy with respect to s are presented in Ap- 
pendix lA II In the vicinity of the criticality, dominant 
singularities of D, n and S are determined by a thermal 
average of </? {{'-p)) from 

D[s\ « d%,*]-^[s\{^) (35a) 
dh 

n[s\ « n[s„*] + |-[sl(^) (35b) 

S[s\ « S%,*] + ^[s\{^). (35c) 

It is remarkable that as elucidated in Eq. H35a|l . the 
relevant scalar field represents a shift of the double oc- 
cupancy D from the critical value at least around the crit- 
icality, since dh/ dU generally does not vanish. Similarly, 
Eqs. Ij35b|) and (|35c|) suggest that the field tp can also rep- 
resents the shift of the electron filling n or the entropy 
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S from the critical value only in the case oi dh/dfj, ^ 
or dh/dT ^ 0, although in the {U, fi, T) space, there ap- 
pears a curve of dh/dfi = 0. 

Except the special curve determined by ft, = dh/dii = 
with r < Tc (open dotted curve in Fig. 13 the double 
occupancy D and the electron filling n show the same 
discontinuities as {cp) at the first-order Mott transition 
boundary h = with r < Tc (shaded surface in Fig. [21. 
Namely, a first-order Mott transition takes place whether 
it is bandwidth-control or filling-control. The first-order 
filling-control transition is nothing but a phase separa- 
tion. The discontinuities again follow the scaling rela- 
tions given by Eqs. H34c|l and (|34el) . Such discontinuity 
of D was obtained in previous worka^iiSiSiiS^, and that of 
n was also obtainediSiSiil. 

Toward the curve h = dh/dfji — 0, the strength of 
the singular component dh/dfj, vanishes. Although the 
strength may have a linear dependence on a distance from 
the curve, it depends on details of a structure of the 
metal-insulator phase boundary. 

At the special curve h = dh/d^i = 0, /x is eventually 
decoupled from the field h conjugate to if and hence n 
is decoupled from (ip) up to a linear order. This means 
that only a first-order bandwidth-control transition takes 
place there but not a filling-control. Then, n is contin- 
uous and a phase separation does not occur. Actually, 
an infinitesimal deviation of fj, around the special curve 
always gives a metallic state on the both sides. 

Lastly, we discuss the entropy S. Since the Mott tran- 
sition surface has a slope against a temperature variation 
as already noted, dh/dT does not necessarily vanish on 
the surface even at the critical end curve, and hence there 
appears a discontinuity of S which scales as {ip). This 
leads to a finite latent heat with varying the tempera- 
ture across the first-order Mott transition surface, as in 
liquid-gas phase transitions. 



with X4 being given in Eq. HA7|) . Therefore, the doublon 
susceptibility xdd diverges toward the criticality as ob- 
tained for the Hubbard modetfiSi, and the critical expo- 
nents are the same as those for x given in Eqs. I|34dl) and 
(|34fl) . Similarly, the charge compressibility Xnn and XDn 
also diverge at the critical end curve (thick solid curve in 
Fig. 121) except the special critical end point (?7cr, Mcr, ^cr) 
determined hy Sr = h = dh/dfi = (open circle in 
Fig. 121). Similar compressibility divergence was obtained 
in infinite dimensions'^. At (C/cr, Acr, T'cr), while xdd di- 
verges as already noted, xon and Xnn become regular, in 
contrast to the other points in the critical end curve. This 
is because dh/dfi and thus the second terms of Eqs. (|36b|) 
and H36c|) vanish there. 

Next, we discuss a scaling behavior of Cu,fj., and 
Ljj. These quantities can also be expressed as 



Lu 



,dh 2 , dr d'^h . . 



dh dh 



dr dr 



dT dU^ dT dU^'^ 



dh dh 



dr dr 



-X 



dTd^i^ dTd^j, 



d^h 

' dUdT 
d^h 
dfidT 



(37a) 
(ip) (37b) 
(^). (37c) 



When dh/dT does not vanish due to a finite slope of 
the Mott transition surface with respect to temperature, 
singularity of C(7,/x is dominated by the first term in 
Eq. (|37a|) and hence the critical properties are determined 
by those of x given in Eq. (|34|) . as in liquid-gas phase 
transitions'^^. Namely, 

C(7,p|r,/i=o |r - TcrT'' oc |T - Tci-T'^ (38a) 

Cu,f.\r^r...k (X |ftri+i/^cx|r-r„ri+i/^(38b) 

Otherwize, it reduces to the conventional scaling rela- 
tions for the specific heat with h being fixed constant; 



D. Second derivatives of free energy 

Here we reveal critical properties of susceptibilities 

rp ( d^F 
J 1 gj^ 



at 



-T 



XDD, XDn and Xrm, Specific heat Cu,n 
constant U and /i, and latent heats 

r(|f),,,andL^ = -T(^)^=T(f^^^ 

sions for general second derivatives of the free energy with 
respect to s are presented in Appendix lA 21 

First, we discuss susceptibilities. Their singularities in 
the vicinity of the criticality are given by 



[dTduJ^ 

Expres- 



XDD 



XDn 



,dh 2 , dr 2 



dh dh 
diidU^ 



dr dr 
dUdJl^^ 



d^h 
dip 

d^h 



,dh.2 .dr.^ 



dUdfi 
dfi' 



(ip), (36a) 
(^), (36b) 
(36c) 



Ch\r.h=o oc Ir-rcrp" oc IT-Tcrp" (39a) 
Ch\r=r^,,h (X \h\-"/^' ^\T-T,,\-^/^'. (39b) 

The same scaling properties applies to and Lu, ex- 
cept at (J7cr, Acr, ^cr) whcrc Lu becomes regular due to 
vanishing singular components. 

It is remarkable that the obtained divergence of spe- 
cific heat Ch at the critical end curve for d < = 4 
contrasts to MF arguments such as the DMFT— and the 
Brinkman-Ricai^ or slave-boson MF— theory. In the MF 
universality which is indeed realized in d > dc, the spe- 
cific heat and thus the quasiparticle effective mass do not 
diverge for T > even at criticality. This reflects that 
the inverse of the quasiparticle renormalization factor. 



l/z(k) = l-^(c. = 0,k). 



(40) 



remains finite at the finite-temperature criticality. As 
we have obtained above, diverges at the finite- 
temperature critical end curve for 2 < d < dc- The 
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diverging quantity is, however, not l/z(k) but the in- 
verse of the suppression factor of the Fermi velocity due 
to spatial fluctuations. 



«'(k)K(k) 



1 



dReY. 



{lu = 0,k)/w(k) 



(41) 



with k± being the component of k perpendicular to the 
Fermi surface. For d = oo, v'^(k)/v^(k) — 1, since there 
exists no k dependence of the self-energy part; a regular 
behavior of Ch in the case of d = cx) is immediately under- 
stood as a typical MF behavior. We can obtain a further 
insight into the critical properties of the Fermi velocity 
z;''(k) = z(k)i;J:.(k) in finite dimensions. The entropy S 
of the quasiparticle system is given by the following k 
integral on the Fermi surface Sp, 



S(xT 



datt 



1 



{2Tr)d^^ z(k)w^(k)' 



(42) 



where ak represents a (d — l)-dimensional infinitesimal 
surface element around k. We assume that both w^(k) 
and Wfe(k) vanish most rapidly toward a certain Fermi 
momentum k^- as Sk'^, where 6k± is the amplitude of 
the tangential component of k — k^ to the Fermi surface. 
Since from Eq. (|39|) a scaling dimension of S is (a—l) jv = 
—d + 1/v and z(k) is finite, then v{k) should have a 
scaling dimension oi \/v + vj — 1. Namely, 



w'Xk) = ^(k)t;I.(k) cx T^^" « |T - 
C l-v{l- vj). 



(43) 
(44) 



Particularly, in the case of a linear dispersion relation 
tTJ = on the whole Fermi surface, toward the criticality 
along the metal-insulator boundary h — 0, v'^{'k) vanishes 
as l/ln|T - T*J in two dimensions and as \T - T*^.\°-^^ 
in three dimensions. If the bare Fermi velocity vanishes 
as 6k\\ as in a van-Hove singularity, namely, w = 1, then 
v'{k) vanishes as Sk\\\T ~ T*J in 2 < d < 4. These 
sharply contrast to the case oi d > dc where v'-'/v^ does 
not vanish. 



E. Discussion 

As already noted*' the same singular behaviors as 
((/?) or X can be observed in many other physical quan- 
tities for transport, magnetic and optical properties. 
However, in spite of intensive and extensive experimen- 
tal investigations on the Mott transition systems like 
V2-.(Cr/Ti),03-2„ NiS2-.Se„ ac-(BEDT-TTF)2X, etc, 
there are only a few systematic experimental studies of 
the thermal criticality of Mott transition^^. Here, we dis- 
cuss possibilities of experimentally verifying the critical 
behaviors from magnetic, transport, optical and spectral 
measurements. 

First, we discuss the dc longitudinal resistivity p. 
Around the Mott criticality, it can be expanded in (ip). 
Then, p and dp/ds follow the same scaling relations as 



(</?) and X, respectively, in other words, D and xdd- This 
is consistent with recent experiments on V2-a:Crj;03 by 
Limelette et alm^. Similar discontinuity of p was also 
found in NiS2 though scaling analysis has not been per- 
formed'*^ . 

We proceed to magnetic properties. In the case of 
single-band Hubbard model, the double occupancy D 
and the electron filling n directly determine the local mo- 
ment (S'foc^)^/^ = (l/2)Vn- - 2D. Therefore, if spin cor- 
relation length is so short that the uniform spin suscep- 
tibility xss follows the Curie- Weiss law, then the Curie 
constant should be proportional to {n — 2D) /A and hence 
directly show the singularities of D and n^. Namely, 
the Curie constant exhibits a discontinuity at the first- 
order transition, whether we approach the first-order 
Mott transition surface in Fig.|21by decreasing the carrier 
doping concentration or changing U /W . More generally, 
Xss can be expanded in terms of {lp). It discontinu- 
ously jumps at the first-order transition and dxss/dT 
diverges at the criticality. The discontinuity was experi- 
mentally observed, though scaling analysis has not been 
performed. The critical exponents are given by those for 
D and xdd- The same critical properties as xss apply 
to the Knight shift and the spin-lattice relaxation rate. 

In order to probe the detailed nature at the first-order 
Mott transition, NMR measurements are useful when 
there exists a sufficiently strong coupling of the nuclear 
spins to the electron near the chemical potential^. By 
observing a splitting of the NMR/NQR spectral peak 
into two, one can detect an emergence of a macroscopic 
inhomogeneity and coexistence of the metallic and the in- 
sulating states at the first-order Mott transition bound- 
ary, which have different values of D, and even n when 
the filling-control Mott transition also occurs. It is more 
likely that in the time scale of the NMR, macroscopic 
coexistence disappears and the globally stable state is 
realized. In this case, the peak discontinuously jumps. 
This may be experimentally verified in the future. So far, 
NMR measurements in V2-a;(Cr/Ti)a;03_j;^ were done 
mainly at low temperatures, but neither around the Mott 
critical end point nor along the first-order Mott transi- 
tion boundary. Alloying for a carrier doping into Mott 
insulator induces disorder effects. Nevertheless, around 
the thermal Mott criticality, the above critical properties 
should dominate over the disorder effects. 

The divergence of the charge compressibility should be 
observed as the dip of the sound velocity via ultrasonic— 
and Brillouin scattering measurements. This reflects that 
the electronic contribution to the sound velocity, which 
is inversely proportional to Xnm vanishes at the critical 
end curve except at (C/cr, Mcr, ?cr)- Fournier et al. has 
observed the critical enhancement of Xnr^ by the ul- 
trasonic measurements, though scaling analyses have not 
been performed. 

Now we discuss spectral properties such as spectral 
weight of coherent one-electron excitations and Drude 
weight. As in d = oo, {ip) also measures the coher- 
ent spectral weight w observed in photoemission; w = 
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Wcr + w' (if) around the criticahty where Wcr is a value of 
w at a critical point and does not vanish in general. This 
property of w can also be detected through optical mea- 
surements. In many cases, the optical conductivity (t{lo) 
nearly splits into low-energy coherent and high-energy 
incoherent spectral weights. While the latter part ap- 
pears above the Hubbard band gap energy, the former 
has a total weight proportional to w'^. In particular, the 
Drude weight for the lowest-energy excitations has also 
the same critical properties as D. 

In the rest of this section, we discuss implications of the 
present theory. The distinction of metal and insulator is 
strict at T = where the dc conductivity can vanish in 
the insulator. On the contrary, at T > 0, a small but 
nonzero conductivity appears due to thermal activation. 
This makes it difficult to identify the insulator at finite 
temperature. Nevertheless, it is meaningful to categorize 
the Mott insulator as a distinct state from metal as in 
the present arguments. This is possible, because around 
integer but partial electron fillings, two distinct states, a 
globally stable and a semi-stable ones, appear near the 
first-order transition terminating at criticality. We have 
shown that these states, namely, a strongly correlated 
Mott insulator and a less correlated metal, correspond 
to the oppositely polarized configurations of the Ising 
model, as discussed by Kotliar et al. for d = cc-~. 

The universal properties obtained for I?, n, xdd, XDu 
and Xnn implies that if one could construct a GL free 
energy expansion in terms of two fields and 0„ rep- 
resenting D and n by appropriately treating the local 
constraints in a slave-particle representation, then the 
coefficients to (p'jy, 4'D4'n and (jj^-^ should vanish similarly 
at the criticality, except the special bandwidth-control 
critical point (J7cr, Acr, ?cr) where the field ipn becomes 
massive. Furthermore, the present theory poses a severe 
constraint on unified analyses of bandwidth-control and 
filling-control Mott transitions; if the bandwidth-control 
Mott transition is of the first order at the open dotted 
curve h = dh/dfj. = in Fig. [3 the Mott transition at 
the shaded surface are also of the first-order, whether 
the control parameter is D or n. This reflects that two 
fields and 0„ appearing in such analyses in terms of 
the slave-particle representation have the same origin of 
the charge degrees of freedom and must be described by 
a scalar field if. Such free energy functional in terms of 
both and is still not easy to handle due to difficulty 
in treating the local constraints. 



V. INTERPLAY BETWEEN THE INHERENT 
MOTT TRANSITION AND OTHER 
SYMMETRY-BREAKING TRANSITIONS 

So far, we have considered a prototype of Mott transi- 
tion occurring without any spontaneous symmetry break- 
ing, which has a phase diagram shown in Fig. ^ (a)- In 
this section, we discuss its interplay with other sponta- 
neously symmetry-breaking transitions. 



A. Mean- field analyses 

We concentrate on the interplay between the scalar 
field ip and the 0(TO)-symmetric vector field (f) represent- 
ing the charge and other degrees of freedom, respectively. 
For instance, m = 1, 2 and 3 describe interplay with 
charge ordering, superconductivity and magnetism, re- 
spectively. Then, according to Eq. 1)25(1 . our GL free 
energy can be written in the form 



GL 



T 



hifii^) + i (r^^(x)2 + (V(^(x))2 



1 / 



+ ^ (r^0(x)2 + (V0(x))2 j + |^(x)0(x)2 



' 4 



(45) 



r^p and are the second-order coupling constants for Lp 
and (/), respectively, and denote the third-order 
coupling constants for (/? and the mixing of (/? and 0, re- 
spectively, u^^ and represent the fourth-order cou- 
pling constants of (/? and </>, respectively, while u^p^ that 
of their mixing. Within the GL theory, phase diagrams 
of the above free energy are given by Fig. ^ (b-d). 

In the case of Fig.n(b), the inherent Mott transition, 
which occurs without any spontaneous symmetry break- 
ing, is removed by emergence of the antiferromagnetic in- 
sulating phase, and only a Mott metal-insulator crossover 
remains at the high-temperature phase. This occurs 
when - r^u^^/u^^ > holds at /i = -gr^/i2u^^) in 
the ordered phase < 0. Physically, this corresponds to 
the case where a partially filled band disappears in the re- 
duced Brillouin zone of this antiferromagnetic phase and 
hence it smears out a phase transition only associated 
with the charge degrees of freedom. Then, this phase 
transition is an Slater-type Mott transition. The phase 
diagram appears in d > 2 when a geometrical frustration 
is so weak that it becomes much larger than the Mott 
temperature scale determined by = 0. 

In the case of Fig.n(cl) or (c2), for instance, antiferro- 
magnetic phase transitions take place on the both sides 
of Mott transition. Antiferromagnetic metal and anti- 
ferromagnetic insulator are separated by the first-order 
Mott transition at low temperatures. This occurs when 

± g\/ -r^ju^^ - r^u^^ju^p^p becomes negative with 
decreasing temperature below the Mott critical temper- 
ature along the first-order Mott transition boundary. 
Discontinuity of {ip) at the first-order Mott transition 
boundary introduces a jump 2gyJ^^7\^Ju^ in the mass 
term of (j), producing a difference in transition tempera- 
tures on the both sides of the first-order Mott transition. 
The phase diagram (cl) has been observed in NiS2-3;Sea; 
and i?Ni03 for interplay with AF, where antiferromag- 
netic metal and insulator are separated by the first-order 
Mott transition at low temperatures. Physically, this cor- 
responds to the case where effects of the geometrical frus- 
tration increase from the case of (b) and the Mott crit- 
icality appears above the Neel temperature. The phase 
diagrams in (c) are also relevant to the case where there 
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appears the interplay of Mott criticality with multicriti- 
cal phenomenon of SC and AF^ as in k-(BEDT-TTF)2X. 

In the case of Fig. ^ (dl) or (d2), an antiferro- 
magnetic phase or a superconducting phase appears 
at low temperatures in the insulating or the metal- 
lic side of the first-order Mott transition. This oc- 
curs when only - \g\^ -r^ju^^p - r^u^^/u^^ but not 

Tcj, + |g|\/-vA<w ~ r^u^^/u^^ turns to be negative 
with decreasing temperature below Tcr along the first- 
order Mott transition boundary. In the issue of inter- 
play with antiferromagnetism, this corresponds to sys- 
tems with geometrical frustration further amplified from 
the case of (c) , smearing out the antiferromagnetic metal. 
The phase diagram (dl) has been observed in such sys- 
tems as \2~xGt:xOj^ and k-(BEDT-TTF)2X under a 
magnetic field^^. The phase diagram (d2) is relevant to 
K-(BEDT-TTF)2Cu2(CN)3^3 for interplay with SC and 
i?2Mo20'^ for that with ferromagnetism. 

Within phenomenological GL arguments, it is also pos- 
sible that a first-order transition and its critical end point 
emerge in the ordered phase, though it is not shown in 
Fig- n This occurs when — rci,u^ci>/u^^ becomes neg- 
ative along the curve of ft. = —gr(f,/U(f,(f, below a certain 
temperature in the symmetry-broken phase < 0. To 
our knowledge, such phase diagram has not been experi- 
mentally observed in correlated electron systems. It may 
occur in charge- ordered spin- disordered phase of quarter- 
filled systems, where once the staggered charge order- 
ing appears then the system effectively corresponds to a 
half-filled one with the spin symmetry being unbroken. 
This phenomenon can also be understood by the present 
framework from the following reason: and <f> represent 
a uniform and a staggered charge densities, respectively. 
While the Hamiltonian does not generally contain a sym- 
metry under ip —ip, it is invariant under (j) — s- —0. This 
leads to the free energy of Eq. 145|) . 



RG equations are obtained as 



B. Fluctuation-induced first-order phase 
transitions 



If the critical regions of two degrees of freedom overlap 
appreciably in the phase diagrams, then interplay be- 
tween these critical fluctuations often alter structures of 
the phase diagrams. Below, taking into account critical 
fluctuations by means of one-loop renormalization group 
(RG) analyses, we show that fluctuation-induced first- 
order transitions as multicritical phenomena may occur 
in the vicinity of the meeting points of the transitions. 

We introduce a momentum cutoff A = e"^, and define 
dimensionless coupling constants h = A^^^'^/^/i/y^, 

= A-^I^r,, g = A-^+d/z^^ ^nd u,^ = A-^+'^Iau,^ 
with Id = 2[(47r)''/2p(^/2)]-i. Then, we perform renor- 
malization group (RG) analyses. Up to one-loop order. 



dg 
dX 


-\- 1777, -|- 2 -|- (777, - 

-g^{l + mup^/{6u^p))]g 


1- 8 )7y,^^ -\- 

1 ^ j Ll/(p(p 1 


± 1 { ^llirs.r^ 

(46a) 


dX 


-eu^^ + 9ul^ -f mul^ + 




(46b) 


du^^ 
dX 


-£U00 + (m + 8)ul^ -\- u\ 




(46c) 


duip^ 
dX 


[-£ + 3m^^ + (m + 2){t0^ 




P0 + 5"^! 



(46d) 



with £3 — u 2 



3-|>0ande = 4- d>0. 
There exist three fixed points of (g, w^^); a 
decoupled Ising-MF fixed point (0, e/9, 0, 0), a decoupled 
Ising-O(m) fixed point (0, e/9, e/(m -I- 8), 0, 0) and the 
0(m + l)-symmetric fixed point (0,e/(m -I- 9),e/{m -\- 
9), e/{m -\- 9)). They are all unstable to an infinitesimal 
value of g. The two decoupled fixed points are character- 
ized by a simultaneous emergence of the inherent Mott 
criticality for ip and 0(m)-symmetric one for (f) which 
are decoupled from each other. It is also unstable to 
Uip4>, Uipip ^ 0. The 0(m -f- l)-symmctric fixed point leads 
to a bicritical behavior and is only partly stable within 
the space of (wi^i^, Wi^^, m^^). We note that up to the 
one-loop order, for m = 2 or 3, there does not exist a 
biconical fixed point which usually leads to a tetracriti- 
cal behavior. In Fig. we show RG flow trajectories of 
the parameters {uipip/u^^, u^(j,/up^, g/u^^) in the case of 
m = 3, d = 3 and Uif,^ > 0. 

In the plane of g'^/ucp4> = 0, the region nearly deter- 



mined by u, 



> 1 shrinks into the 0{m -f- 1) 



symmetric fixed point as RG fiows show. The other 
part of the plane flows toward Uipip/utp^ —00 and/or 
U(f>cf>/u^ip —00. Such runaway flows toward negative 
values of u^p^p and/or indicate a fluctuation- induced 
first-order phase transition^. Usually, a tricritical be- 
havior appears in the case of runaway flows. 

A flnite third-order coupling g renders the interplay 
between p and 4> more severe, which is also displayed 
in Fig O The flows are largely categorized into three 
depending on the bare coupling constants: 

(I) If the bare coupling constants satisfy 
u^uyUsslu^pS ^ li the parameters {upp,Up^,u^^) 



first flow toward the 0(m -|- l)-symmetric fixed point 
and then quickly turn to go to infinity. Then, 
becomes negative and sign changes of Uipp and u^^ 
follow, although it is not shown in the figure. In this 
case, the system undergoes a fiuctuation-induced first- 
order phase transition to a phase which is likely to have 
both long-range orders {p),{<j)) ^ 0. However, detailed 
properties of the ordered phase actually depends on 
higher-order terms which are not included in Eq. (|45|l . 

(II) If the bare coupling constants satisfy 
UipipU(t,(t,/v?p^ <C 1, the parameters flow toward in- 
finity. In this case, Uipp and/or turn to be negative 
while u^(j, remains positive. This leads to a fluctuation- 
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induced first-order transition to a phase which is likely 
to have only an long-range order of either Lp or (p. 

(Ill) We consider the cases where the bare coupling 
constants satisfy UipipU^^/u^^ ~ 1. In this case, the pa- 
rameters (uiptp,u^^,Ucj,cj,) first approach the 0{m -|- 1)- 
symmetric fixed point, and the flow is modified by the 
fixed point into a nearly vertical one. Then, while the 
relation Uj^j^ ~ u^^ ^ u^^ is approximately satisfied, 
g'^/u^pci, rapidly increases. This gives rise to a bicritical 
scaling behavior away from the transition. Finally, the 
parameters fiow to infinity and the fiuctuation-induced 
first-order phase transition noted in either (I) or (II) ap- 
pears instead of the bicriticality. 




FIG. 3: (color online) RG flow diagram in the space of 
Uti,4,/u^4,, g /uifi4,) for m = 3 and d = 3. Dashed 
curves represent RG flows in the plane of /u^p^ = 0. They 
are separated by the two dotted curves. Solid curves represent 
RG flows for (j^ juip^ > 0. 



VI. SUMMARY 



the chemical potential ^ and the temperature T, there 
appears a first-order Mott transition surface terminating 
at a critical end curve. Around the criticality, the dou- 
ble occupancy, the electron filling and the entropy have 
(y5-linear terms. They show discontinuities and diverging 
responses at the first-order and the second-order Mott 
transitions, respectively. The divergence of specific heat 
indicates that of quasiparticle effective mass at the criti- 
cality. There exists a curve where the Mott gap closes. At 
this curve, a filling-control transition does not occur but 
a bandwidth-control, the (/s-linear term in the electron 
filling vanishes, and the related singularities disappear. 
We have also provided principles for experimental tests 
of the proposed critical scaling behaviors and emergence 
of a phase separation from transport, magnetic and spec- 
tral measurements. 

Various experimentally observed phase diagrams have 
been discussed in terms of the GL theory for interplay be- 
tween the inherent Mott transition occurring without any 
symmetry breaking and a symmetry-breaking phase tran- 
sition such as charge ordering, superconducting and an- 
tiferromagnetic transitions. It has also been shown that 
Mott transition accompanied by a symmetry breaking is 
described by the symmetry-breaking phase transition un- 
der renormalization effects in the proximity of the Mott 
criticality. In the presence of two critical degrees of free- 
dom, there is a possibility of observing a bicritical phe- 
nomena as a crossover. If the renormalization becomes 
strong, even a fiuctuation-induced first-order transition 
occurs as a multicritical phenomenon. 

Model calculations of the obtained critical properties 
and an application to quantum Mott transition are left 
for future studies. Here, we only note that present ar- 
guments can not be directly extended for addressing is- 
sues of quantum Mott transition and quantum-classical 
crossover. This is because at T = 0, susceptibilities 
Xdd, XDn and Xnn vanish in the Mott insulating states, 
in contrast to the case of T > 0. The present frame- 
work is also applicable to other correlation-driven metal- 
insulator transitions, like those associated with a charge 
ordering phenomenon. 



We have developed Ginzburg-Landau theory of classi- 
cal Mott transition for the dimensionality d < oo by reex- 
amining the framework of the cluster dynamical mean- 
field theory. Then, we have confirmed that the ther- 
mal critical phenomena for the Mott transition occur- 
ring without any spontaneous symmetry breaking can be 
described in terms of a liquid-gas phase transition char- 
acterized by the Ising universality class, in accordance 
with a recent experimental scaling analysis of the resis- 
tivityS^. There appears a relevant scalar field ip repre- 
senting the charge degrees of freedom characterized by 
the double occupancy. This plays a role of "order pa- 
rameter" and describes a bosonic soft mode at Mott crit- 
icality. In the space of the local Coulomb repulsion U, 
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APPENDIX A: U, n AND T DERIVATIVES OF 
THE FREE ENERGY 



1. First derivatives 

Here, we give an expression of the first derivative with 
respect to Sj which is a jth component of s; 



I dF 1 

AT a„ l-'^-l AT" 



dF, 



GL ( 



Nc dsj 

nJ 

1 dF, 



Na' dsi 



dh 



1 



1 du 



ds 



2ds 



(2)r 



Neds 



(Al) 



where 



X3 

vanishes at criticahty and 



1 d^F 
Ndhdr 



1 d^F 



X4 



N dr^ 

has a critical divergence characterized by 



X4|/i=0,r OC 
X4:\h,r=rar 



\h\ 



(A6) 



(A7) 



(A8a) 
(A8b) 



where (• • ■) denotes a thermal average with respect to 
-Fgl, and in the case of temperature derivatives, we have 
omitted some sub-leading terms. In the case of h deriva- 
tives of the free energy, h and 5r dependences appear in 
(tp) but not in their expansion coefficients. In the present 
case of Sj derivatives, however, the expansion parameters 

(i) . 

alj with i = 1, 2, • • • depend on the original model pa- 
rameters s as 



with the specific heat exponent a — 2 
subdominant for a < 7. 



dv and hence is 



APPENDIX B: DERIVATION OF 
RENORMALIZATION GROUP EQUATIONS 



.(1) - 



.(2) - 



dh 
dsj 
1 dr 
2d7j 



du 1 
'ds'jll, 



4-3 



du 1 
'ds'jWc 



dx(Kx)-(^))3) (A2) 
dx(Mx)-(^))2)(A3) 



They satisfy daf-' /dsj' — da^9 /dsj, which is referred to 

j 

as b'-'^K In particular, at the Mott transition boundary or 
a metal-insulator crossover h = 0, we obtain 



,(1) - 



^(1) = _ 



dh 
dsj 

d^h 
dsjdsji 



(A4a) 
(A4b) 



2. Second derivatives 

Here we give an expression of the second derivative 
with respect to Sj and Sj'] 

1 d^F 



Nc dsjdsj 



^ 1 f(9FGL 

TVc ds, 



idFcL^^^dFoL 



.dF 



GL 



dsi 



dh dh ^ dh dr 

dsj dsj' dsj dsj' 



dsj, 
dh dr 
dsj' dsj 



dsi 



» 



d^F 



GL 



)X3 



dsjdsj, 
dr dr 
dsj dsj' 



Here, we derive the renormalization group equations 
given in Eq. H46|l . First, we add to the GL free energy 
Eq. 145|) a third-order coupling term {gip/Q) J dxip{x)^ 
and define the corresponding dimensionless coupling con- 
stant = ^/Tdgip. Then, the RG equations are obtained 
up to the one- loop order as follows; 



dX 
dg_ 

dX 



-£3 + {m + 2)u^,p + (to -I- : 



mf]g (Bla) 



du,, 



dX 



mu^^ - 



1 



(3: 



;;2 



mg'^) 



(Bib) 
(Blc) 



dX 



-eu^j,,), + (to -I- 8) 



(Bid) 



du 



v4> 



dX 



V0 + [Sutpi^ + (to + 2)u00 + 4:U^p]u^^ 



+ ~g^ + ~gvf + ~9l~9'' ■ 



(Ble) 



Next, we impose = at each value of A 

^ by renormalizing (p{x) = ^p{x)K^'^~'^^/'^ as (p{x) — > 

g(<i-2)(5A/2^|'2,^ — Sgtp/{6uip^). This introduces an addi- 
f • • -tional term Uip^/{6uipip) ■ dg^/dX to the RG Eq. (|Blbp . 

Finally, substituting = into the RG equations, we 
(A5)obtain Eqs. gH). 
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